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Abstract 



O A set of coupled kinetic equations describing in the Abelian approximation a mix- 

ture of quarks and self-interacting gluons is formulated and solved numerically. The 
O ■ model includes the Schwinger-like mechanism for particle creation in a strong field 

O ■ as well as two-particle elastic collisions between all mixture components in the Lan- 

t^- ! dau approximation of small-angle scattering. The process of equilibration at the 

initial energy density exhibits a dominant quark creation in the very early time of 
; interaction. It is shown that damping of energy density oscillations due to elastic 

scattering of perturbative quarks and gluons is not strong enough to reach thermo- 
q . dynamic equilibrium in a reasonable relaxation time. A possible account for a such 

behavior is discussed. 
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X\ ■ 1 Effective Lagrangian 



Quark-gluon plasma (QGP) evolution is considered in a following physics scenario. In early 



stage of heavy-ion collisions at ultra-relativistic energy a chromo-electric flux tube may be 
stretched between nuclear residuals and this strong field results in a spontaneous vacuum 
pair creation (Schwinger mechanism). In its turn, the created particles can influence on this 
background field (back reaction process) and also suffer subsequent rescatterings relaxing 
to some equilibrium state, in general. This picture is treated in the Abelian approximation 
based on the kinetic equations (KE) derived in [p], |2| . In our paper we focus on the relative 
role of components of quark-gluon mixture in the pre-equilibrium stage as well as on the 
question how fast the relaxation process is. 

To describe QGP we use the following effective model Lagrangian ||: 

L(x) = d^*(x)d^(x) -m 2 + |$(x)| 2 - e\ |$(x)| 4 
+te + {[V n $*(x)]$ 2 (x) - [$*(x)] 2 V n $(x)} 

+ -?/>(x)7 M Of, ip(x) - m-ip(x)ip(x) 

+ y ($(x) + $*(x)) ij(x){-fn)i/j(x) , (1) 

where indexes '+' and '— ' are used for bosons (gluons) and fermions (quarks), respectively, 
and V„ = n M (9 M . 
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We suppose that the total bosonic field can be decomposed into a mean-field contribution 
$>o(t) and its fluctuations 

$(x) = $ {t) + <p(x). (2) 

We consider <3> as a neutral space homogeneous background field. The field of fluctuations 
<p(x) is complex and corresponds to charged field with vanishing mean value. The 4 term 
in our Lagrangian (|l|) simulates the gluon self-interaction. Keeping only the second order 
terms in fluctuations, the corresponding equations of motion read 

[iyfy - e_7%$ (t) - m_] ip(x) = 0, (3) 

<p{x) = 0, (4) 
d 2 Mt) + rnl$ (t)+4e 2 + ®l(t)-j-(t)-j + (t) = 0, (5) 

where the covariant derivative = + ie + n M $o(^)- 

Alongside with the Dirac equation (§) and the Klein-Gordon equation @ we get eq.© 
for evolution of the background mean field. One sees that equations of motion are self-consis- 
tently coupled: the created particles generate the currents j±(t) which form the background 
field according to ([5]). In the mean- field approximation these currents are given as 

3-if) = -e_ < ^0X7^0) >> (6) 
j+{t) = -e+ < iip*{x)V n (p{x) - i(V n (p*{x))(p(x) 

-2e + %(t)cp*(xMx) > . (7) 



d;d» + m 2 + 



2 Kinetic equations 

Starting from these equations of motion one can introduce quasi-particles and obtain the KE 
for the single particle distribution function on the basis of diagonalization of Hamiltonian 
by the time-dependent Bogoliubov transformation. This result is exact in the mean-field 
approximation for a space homogeneous field. In our model, the vector n M corresponds to 
Flux Tube Model (FTM) geometry and is chosen to be (0,0,0,1). So, we arrive at the 
following KE: 



df±{p,t) f ^ df±(p,t) Q n 

+ e±a{t) — q = S±(p, t) + C±(p, t) . 



Of 



dp\\ 



(8) 



Here, the chromo-electric field strength is a(t) = —d<&o/dt, C±(p,t) is the collision integral 
and S±(p,t) is the source term: 



i r 

S±(p,t) = -W ± (p,t) J dt'W ± (p,t,t') x 



[l±2/ ± (P±(*,0,0]cos |2 / dTu±(p;t,r) 



(9) 



with the transition amplitude 



W ± (p;t,t') 



e±a(t')P ± (t,t') 

<4(p;*,*0 



p ± (t,t>) 



-1 2s 



(10) 
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where s + = for gluons and s_ = 1/2 for quarks. The quasi-particle energy entering this 
equation is 



ul(p,t,t') = e 2 ±x + Pl{t^ 



where P±(t, t') = p\\ — e± J dra{r) and e\ ± = m\ + p\. General features of this source term 

v 

are described in || £|. 

Kinetic equation (§) may be rewritten in the form of a system of ordinary differential 
equations [|, [|: 



df±(p,t) 
dt 

dv±(p,t) 
dt 

du±(p,t) 
dt 



-W±(p,t)v ± (p,t) + C±(p,t), 



W±(p,t) (1 ± 2f±(p,t)) - 2uj ± (p) u±(p,t), 
2u±(p) v±(p,t), 



where two new functions have been introduced : 

/ 

Jo 



(12) 
(13) 
(14) 

(15) 
(16) 



u±(p,t) = J o dt'W ± (p,t,t')(l±2f ± (P ± (t,t'),t'))sm[2j ti dru ± (t,T)], 
v±(p,t) = ^dt'W ± (p,t,t , )(l±2f ± (P ± (t,t , ),t'))cos[2^dru ± (t,T)}. 

3 Collision integral 

Neglecting particle-particle collisions does not result in complete picture of a collision at 
large time moments when the interaction force between particles is getting greater than 
the mean-field interaction. However, the direct evaluation of the CI C±(p,t) gives rise to 
a huge numerical problem. In earlier papers the CI was introduced in the relaxation time 
approximation with time- and momentum-independent relaxation time what does not allow 
to restore true dynamics of the relaxation process and to estimate properly the relaxation 
time 0, Here, the CI will be obtained on a dynamical basis but in a simplified manner by 
making use of the Landau approximation, i.e. assuming small momentum transfer in elastic 
qq-, qg- and gg-collisions. The appropriate cross sections are calculated in the perturbative 
approximation 0: 



da, 



dt 
da 



2 



'.!'.! 



dt 
da, 



Air 



N 2 

c 



s 2 + u 2 



t 2 



s 2 + t 2 



u 1 



2f_ 

3tu 



N 2 - 1 s 2 



ut us st 



t 2 



u- 



'/.'/ 



dt 



2k4(s 2 + u 2 ) 
s z 



1 



N 2 . 



2 N 2 su 



(17) 
(18) 
(19) 



where a s is the QCD coupling constant and s, t, u are Mandelstam's variables. 

As is seen, the cross sections are divergent at small momentum transfer. So, only leading 
terms of the 1/t 2 order are kept in (|IT|)-(|T9"|) which dominate for gluon/quark collisions 
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within the Landau CI approximation. In this approximation one can obtain the Landau- 
like CI: 



C a {fa, t) 







dp aa b 



L -f b (l±f b )-P L fa(l±fa) 



PPa(3 OPba 

where B is a kernel. This kernel is further simplified for massless particles 

B af 3(p a ,Pb) = £ab [(1 - V a V b )S af3 + V aa V bf3 + V ba V af3 } , 



(20) 



(21) 



where v aa = p aa /u(p a ,t), £ ab = 27ia 2 G ab L. Coefficients G ab are defined by corresponding 
cross sections: 



G 



Ml 



G 



<l!l 



G 



8iV c 2 • 



(22) 



4 Calculation results 

The system of three differential equations (|T2])-(|T1|), describing self-consistently the evolu- 
tion of distribution functions and fields, is closed and can be solved numerically as the initial 
Cauchy problem. We choose the zero initial conditions for distribution functions of bosons 
and fermions and for a non-zero initial value of the chromo-electric field strength eait = 0). 
These conditions correspond to the FTM, where colliding ions generate a large-amplitude 
field at passing ions through each other. This field is decaying due to creation of quarks 
and gluons. 

The analysis of abundance evolution within the pQCD model shows the dominant role 
of gluons [0. But as is seen from Fig. la, it is not the case for the model considered : 
Quark production dominates over gluons during few first fm/c in the very early stage of 
interaction. This "fermion dominance" originates from spin effect in the source term ([J) 
suppressing low P\\ gluons |]]. In the early stage of a collision when the created particle 
density is still small, in the phase space there are many free states for the final state of 
created particles and the influence of the Pauli blocking is not so essential. However, with 
subsequent density increase, the fermion creation is suppressed by the factor (1 — 2/1) due 
to occupied states and bosons are enhanced by the (1 + 2/ + ) factor. Therefore, only at 
later time the system evolution resembles that in the pQCD model exhibiting the dominant 
gluon production. 

The mass ratio m + / m_ dependence of the fermion dominance time 7/ is of great interest. 
As shown in Fig. lb, 7/ ~ const in the range of m + /m_ < 0.5 corresponding to non-zero 
fermion dominance time even in the limit of zero boson mass. However, r/ — ► oo when the 
boson mass equals to or exceeds the fermion mass. 

Unfortunately, the CI (]2ll|), taking into account mainly hard partons, is of minor im- 
portance and can not result effectively in quantum oscillation damping of the mean field. 
The mean field causes a rippling excitation of the distribution function and makes impos- 
sible to achieve equilibrium in QGP in the small momentum-transfer approximation used. 
Fig. 2 illustrates this fact for pure gluon plasma and for QGP. In the first case considered 
the system is close to equilibrium, but taking into account additionally the quarks degrees 
of freedom, we obtain far-of-equilibrium dynamics. The latest work shows that this 
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Figure 1: Time dependence of number density for fermions (solid line) and for bosons 
(dotted line) (a) and the fermion dominance time r/ as a function of the m + /m_ ratio (b). 
All results are presented for the system with the initial field strength ea(t = 0) = 10 GeV A . 



problem is not conditioned by using the Landau CI and can not be solved by including 
higher terms in evaluation of corresponding cross-sections. Due to the fermion dominance 
at the early stage, quarks cause significant changes in the gluon distribution function and 
the system evolves extremely slowly towards a quasi-equilibrium state. 



P, GeV 4 P GeV 4 




t, fm/c 

(a) (b) 



Figure 2: Time dependence of longitudinal (solid line) and transverse (dotted line) pressure 
for pure gluon system (a) and that for different parton components in quark-gluon plasma 
(b). 



5 Conclusions 

The self-consistent system of equations for describing early stage of heavy ions collisions 
has been derived. Being solved numerically, the set of equations does not result in a 
fast attainment of a quasi-equilibrium state for the system under discussion. To get a 
reasonable estimate for the relaxation time of quark-gluon plasma one needs to increase the 
cross sections in few times what effectively would correspond to accounting for radiative 
pertubative processes and scatterings of nonperturbative partons, as well. Instead of the 
gluons dominance predicted by pQCD models, we observe the quark dominance in the very 



5 



early stage of particle production at the time scale of few fm/c corresponding to QGP 
evolution. If expansion of excited matter is included into consideration, this time Tf should 
be even longer. 

One should stress that the KE method used is rather simplified with respect to the CI 
( 20|) in KE (§). This Boltzmann-like integral is obtained in the nearest order approximation 
of the gradient expansion for the distribution function, neglecting all coherence effect to 
be related to the mean field. As expect, the consecutive dynamical approach to the CI 
problem can lead to an important correction of the result discussed : Strong quasiclassical 
fields have to influence on quantum fluctuations of quark and gluon fields. This problem is 
a real challenge to the kinetic theory of a particle-antiparticle plasma in a strong field. 
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